arXiv:astro-ph/0512055vl 2 Dec 2005 


Probing for evolutionary links between local ULIRGs and QSOs 

using NIR spectroscopy 

K. M. Dasyra, L. J. Tacconi, R. I. Davies, R. Genzel, D. Lutz 

Max-Planck-Institut fur extraterrestrische Physik, 85741 Garching, Germany 

T. Naab 

Univesitatssternwarte, Scheinerstr. 1, 81679, Miinchen, Germany 

D. B. Sanders 

Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 

S. Veilleux 

Department of Astronomy, University of Maryland, College Park, MD 20742, USA 

A. J. Baker 

Jansky Fellow, National Radio Astronomy Observatory 
Department of Astronomy, University of Maryland, College Park, MD 20742, USA 


Abstract 

We present a study of the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs), merging galaxies of 
infrared luminosity > 10 12 Lq. During our Very Large Telescope large program, we have obtained ISAAC near- 
infrared, high-resolution spectra of 54 ULIRGs (at several merger phases) and 12 local Palomar-Green QSOs, to 
investigate whether ULIRGs go through a QSO phase during their evolution. One possible evolutionary scenario 
is that after nuclear coalescence, the black hole radiates close to Eddington to produce QSO luminosities. The 
mean stellar velocity dispersion that we measure from our spectra is similar (~160 km/s) for 30 post-coalescence 
ULIRGs and 7 IR-bright QSOs. The black holes in both populations have masses of order 10'-10 s Mq (calculated 
from the relation to the host dispersion) and accrete at rates > 0.5 Eddington. Placing ULIRGs and IR-bright 
QSOs on the fundamental plane of early-type galaxies shows that they are located on a similar region (that of 
moderate-mass ellipticals), in contrast to giant ellipticals and radio-loud QSOs. While this preliminary comparison 
of the ULIRG and QSO host kinematical properties indicates that (some) ULIRGs may undergo a QSO phase in 
their evolutionary history before they settle down as ellipticals, further data on non-IR excess QSOs are necessary 
to test this scenario. 


Key words: galaxies: formation, kinematics, & dynamics, infrared: galaxies 


Preprint submitted to Elsevier Science 


5 February 2008 




1. Introduction 


In the local Universe, the best laboratories 
for studying violent merging events (considered 
key-mechanisms in driving galaxy evolution) are 
the ultraluminous infrared galaxies (ULIRGs). 
Several studies indicate that ULIRGs trans¬ 
form gas-rich disks into moderate mass ellipti¬ 
cals thr ough merger-in duced dissipative collapse 
(iKormendv fc Sanders!. il992 ; iMihos fc Hernouisti . 
H99ar Photometric analysis of u lt ralum inous 
merger remn a nts by IVeillenx et a, 1.1 ||20O2) and 
IScoville et alJ < 200(11 indicates that most of 
the sources are are well-fit by an elliptical- 
like r 1//4 light prof ile. Specrtoscopic analys is by 
iGenzel et al . d200lh and lTacconi et al ■ (M) indi¬ 
cates that their stellar kinematics resemble those 
of dispersion-supported systems. 

While the end products of galactic mergers 
can be considered understood, several details of 
the merging process (often related to the physics 
of the gas) are still very uncertain, even in the 
local Universe. A pleth or a of n umerical mod 
els iMihos fe Hernauist. 1996; Springel et al 


2005|) and observ ations llSanders fc Mirabel Il996: 


Kim et all 2002 ) indicates that major starburst 


episodes occur after the first galactic encounter 
and can be present after nuclear coalescence, be- 
fore complete relaxat ion sets in. According to 
iDi Matteo et alJ (2005), the main accretion event 
happens after the coalescence. However, the rel¬ 
ative contributions of starburst and AGN to the 
infrared (IR) emission during each phase of a gas- 
rich merger are still unclear. Whether ULIRGs 
may go through a QSO phase after the nuclear 
coalescence also needs to be confirmed. 

One way to investigate the physical details 
of mergers is to determine the evolution of the 
ULIRG kinematic properties as a function of time. 
Hence, we obtained high-resolution, H- and K- 
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band ISAAC spectroscopic data for 54 ULIRGs 
(spanning a wide range of merger phase) during 
our European Southern Observatory Very Large 
Telescope (VLT) large program (ID 171.B-0442). 
Of these sources (at 0.018 < z < 0.268), 23 are 
binary (not fully merged) ULIRGs and 30 are (re¬ 
laxed) remnants. We extract the stellar dispersion 
a and rotational V TO t velocity from the CO rovi- 
brational bandheads that appear in our spectra 
using th e Four ier correlation quotient technique 
feendeii. 1220'); this method provides the line- 
of-sight (LOS) broadening function with which a 
stellar template has to be convolved to produce 
the observed spectrum. 


2. Dynamics of ultraluminous infrared 
mergers 

2.1. Pre-coalescence phase: the conditions that 
trigger ultraluminous emission 


We derive the kinematics of each individual com¬ 
ponent of the binary systems, to investigate the 
conditions under which ultraluminous infrared ac¬ 
tivity is triggered. From the stellar dispersion and 
rotational velocity, we compute the masses of the 
merging galaxies. We find that that ultraluminous 
luminosities are are mainly generated by almost 
equal-mass mergers; the average mass ratio of the 
binary ULIRGs is 1 .5:1 and 68% of the se sources 
are 1:1 encounters 1 Dasvra et al .. fiooi . Statisti¬ 
cally, our result is in agreement with lu minosity ra- 
tios i nferred from imaging studies (e.g. lKim et all 
EiJ. In individual cases though, luminosity does 
not always accurately trace mass due to stellar 
population and extinction effects. 

To investigate whether dynamical heating 
(mainly of the smaller companion ) can s kew our re¬ 
sults, we re-ran the simulations of Naab V Unrkertl 
( 2003 ) with an V-body-plus-SPH code after re¬ 
placing 10% of the stellar mass with isothermal 
(10 4 K) gas. Following the stellar dispersion as a 
function of time, we found that 3:1 mergers may 
appear to be 2:1 at nuclear separations as small as 
those we observe (7.2 kpc on average). Since only 
a minority of our sources belong to these merger 
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categories, our conclusions are stable against dy¬ 
namical heating. Mergers of mass ratio >4:1 ap¬ 
pear not to drive enough gas to the center of the 
merger to generate ultraluminous IR emission. 


2.2. Post-coalescence phase: the evolution of 
the host properties and the end-products of 
ultraluminous mergers 

By studying the kinematics of ULIRGs at 
merger phases that follow nuclear coalescence 
(remnants) and comparing them to the kinemat¬ 
ics of the binary ULIRGs, we find that the mean 
stellar dispersion increases during the coales¬ 
cence phase by 15 kms -1 (namely from 142 to 
157 kms -1 ). Although this increase is significant 
(within the 0.05 significance level of our Monte 
Carlo simulations), it constitutes only a lower 
limit on the dynamical heating that the merging 
galaxies undergo. This is primarily due to the fact 
that the ultraluminous phases are s hort compar ed 
to the total merger timescales (e.g. iMihostl1999 b 

The ratio between rotational and dispersion ve¬ 
locity of the ULIRG remnants, which equals 0.31 
and increases to 0.62 when inclination effects are 
taken into account, indicates that the remmants 
are dispersion-supported systems, resembling el¬ 
lipticals (Es). To investigate what t ype of ellipti¬ 
cals ul traluminous mergers will form. lGenzel et all 
(l2D0lh suggested that ULIRGs need to be placed 
on t he fundamental plane (FP ) of early-type galax¬ 
ies dDiorgovski fc Pavf3 . fl98 7l. In Fig. 1 we present 
the effective radius- host dispersion (i? e ff — <r) pro¬ 
jection of the plane, constructed from our data for 
54 ULIRGs (triangles). Giant (boxy-isophotal pro¬ 
file) Es are shown in boxes, moderate-mass (disky- 
isophotal profile) Es in filled circles, and further 
cluster Es in open circles (see caption for origi¬ 
nal references). Some luminous infrared galaxies 
(LIRGs; 10 n -LQ < L < 1O 12 L 0 ) are presented 
in diamonds. The fact that ULIRGs clearly popu¬ 
late the moderate-mass ( 10 11 Mq) Es part of the 
FP suggests that these two populations are linked, 
while giant Es probably have a different formation 
history. 

From the dispersion velocity of the stars in 
the bulge of our ULIRG remnants, we compute 


their BH mass Mbh using th e Mbh ~ a rela¬ 
tion (e.g. iTremaine et mi20021. To test whether 
the Mbh — er relation can be applied at the late 
ultraluminous phases of a merger, we ran our sim¬ 
ulations (presented in Sect. 2.1) and followed the 
host dispersion and the gas inflowing to the center 
of the simulation as a function of time. We found 
that by the time the nuclei coalesce, these quan¬ 
tities already scale linearly as long as the amount 
of gas that accretes onto the black hole from its 
surroundings (center of t he simulation) is r oughl y 
constant with time (e.g. iDi Matteo et al 1 EqqI). 
The application of the Mbh — & relation to the 
remnants yields black hole masses of the order 
10 7 - 10 8 M q . 

The Eddington efficiency ?7Edd (the ratio be¬ 
tween the Eddington and the dynamical Mbh) of 
the ULIRG remnants is calculated by assigning 
50% of the IR luminosit y to t he A GN. This statis¬ 
tical assumption follows lGenzel et all : 191)sjh who 
found that some ULIRGs are largely starburst- 
while others are AGN- powered, and may cause 
a few individual sources to accrete at super- 
Eddington rates. If we assign 77Edd=l to these 
sources, the mean Eddington efficiency is 0.5. Such 
high 7/Edd values may be an observat ional confirma¬ 
tion o f the theoretical predictions o f fSnriugel et a, 1.1 
ll2005l) and lPi Matteo et al.l d2005lh These authors 
suggest that after nuclear coalescence, the gas in¬ 
fall to the center of the system is so high that the 
AGN may accrete at near-Eddington rates. 


3. Dynamics of Palomar-Green QSO hosts 
and their relation to mergers 

An evolutionary scenario for the late phases 
of gas-ri ch me rgers, originally based on that of 
ISanders et alJ (jl98«b is as follows: after coales¬ 
cence, the IR emission arising from the nuclear 
starburst and AGN-surrounding dust is strong 
enough to reach QSO-like luminosities. However, 
as the dust and gas start clearing out from the 
nuclear region due to AGN winds and supernova 
feedback, the system goes through a short (up 
to 10 s yrs) optically bright phase, before further 
accretion and star formation are finally prevented 
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1 Springel et al l E55a . 

To test this scenario we obtained spectro¬ 
scopic (VLT ISAAC-PI Tacconi, SPITZER IRS- 
PI Veilleux) and imaging (HST NICMOS-PI 
Veilleux) data for a (small) sample o f Palomar- 
Green QSOs llSchmidt fc Greeii j 198^1 as part of 
our QUEST (QSO/ULIRG evolutionary study) 
project. The VLT sample consists of 12 sources, 
most of which are IR-bright (ratio of integrated IR 
to big blue bump luminosity > 0.46; Surace et al., 
2001). If the above evolutionary scenario is realis¬ 
tic, the IR-bright sources are transitional objects 
between ULIRGs and (optically-selected) QSOs. 

The HST NIC MOS imaging a nalysi s of 7 IR- 
bright QSOs by I Veilleux et alJ <j2006(l indicates 
that their hosts have similar H-band host mag¬ 
nit u des to those of ULIRGs. In a similar vein, 
iGanalizo fc Stockton (feoQlil found that the (op¬ 
tical) spectra of most of their IR-bright QSOs 
(selected from the IRAS color-color diagrams) 
are well fitted by an old population and a recent 
(< 3 x 10 8 yrs) strong starburst component. 

From our high-resolution, large-collecting area 
data (obtained under excellent seeing conditions < 
0.5), we succeeded in extracting host dispersions of 
local QSOs from NIR CO bandheads. The average 
dispersion of the IR-bright QSOs is (so far) sim¬ 
ilar to that of ULIRG remnants (~ 160 kms -1 ). 
Thus, their black hole and host masses are also 
of comparable size (10 7 -10 8 and 10 lo -10 11 M 0 re¬ 
spectively). Since IR-bright QSOs have similar lu¬ 
minosities (> 10 12 Lq) to ULIRGs, their accre¬ 
tion rates are also high (~ 0.5 Eddington). For the 
few optically selected sources observed so far, the 
scatter in the dynamical properties is higher. 

In one of the 12 QSOs targeted by our large 
program, PG1426+015, two nuclei separated by 
^4kpc appear in the acquisition image (also seen 
in HST NICMOS data; Schade et al., 2000). Our 
NIR spectroscopy indicates that the nuclei are at 
the same redshift. The optical spectrum of the 
bright comp o nent has very broad emission lines 
llKaspi et all l2000lh verifying that the emission 
originating from this component can be attributed 
to an AGN. Calculation of the IR luminosity of 
the system (from IRAS fluxes) indicates that it 
is a LIRG. The spectra of the QSO and the sec¬ 
ondary nucleus are presented in Fig. 2. The second 


nucleus does not show any indications of strong 
AGN continuum and resembles those of the bi¬ 
nary ULIRGs; therefore we believe that it is star- 
burst dominated. This system constitutes a good 
example of how mergers can simultaneously trig¬ 
ger strong accretion and starburst events. These 
findings contradict the aforementioned scenario 
in that the QSO phase is reached already before 
nuclear coalescence. 

The IR-bright QSOs seem to differ from local 
QSOs that host even more massive black holes, 
such as the 0.1 < z < 0.25 radio-loud (RL) sources 
or their radio-quiet o ptical counterparts (RQC) 
of iDunloo et alJ ll2003l) . The RL/RQC QSOs have 
black hole masses of order 10 8 and 10 9 M 0 that 
accrete on average at rates ~ 0.05 Eddington, and 
that are located in 5 times more massive and (K- 
band) luminous hosts. The positions of IR-bright 
and RL/RQC QSOs on the fundamental plane can 
be clearly distinguished (see Fig. 1); like giant Es 
vis-a-vis ULIRGs, the RL/RQC QSOs probably 
have a different formation history from the IR- 
bright population. 

Our preliminary results seem to indicate evo¬ 
lutionary links between ULIRGs and IR-bright 
QSOs; however, the fact that some IR QSOs 
have prominent spiral hosts (Surace et al., 2001; 
Veilleux et al., 2006) implies that they may have 
a minor-merger origin (different from that of 
ULIRGs in terms of the mass ratio of the merging 
galaxies). Furthermore, IR-bright QSOs are only 
a small and not necessarily representative part of 
the PG poopulation. We need to enrich our sam¬ 
ple with optically selected QSOs in order to derive 
conclusions about the optically bright phase of 
this scenario. 


4. Conclusions 

From the high-resolution NIR spectrocopic data 
we have obtained for local ULIRGs and QSOs, we 
have found that 

- For IR luminosity > 10 12 L 0 to be triggered 
during a gas-rich merger, encounters of compa¬ 
rable mass galaxies are (typically) required. 

- Evolution (increase) of the host dispersion is ob- 
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served as the merger advances from pre- to post¬ 
coalescence. 

- The merger remnants resemble moderate-mass 
(~10 lo -10 n M 0 ) ellipticals. The black holes 
they host are of the order 10'-10 8 M 0 and, on 
average, accrete at high Eddington rates (>0.5). 

- The IR bright QSO dispersions and black hole 
masses, being of the order 10 7 -10 8 Mq, resemble 
those of ULIRG remnants and indicate a possible 
link between the two populations. Our IR-bright 
sources differ from QSOs that host supermassive 
black holes and accrete at low rates. 

- Imaging and spectroscopy of PG 1426+015 show 
that it is a binary system of nuclear separa¬ 
tion ~4 kpc. Already at this early merger phase, 
one of the components is a powerful QSO (with 
strong dust continuum). 
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Fig. 1. The fundamental plane of early-type g alaxies — c r proje ction). T he giant boxy a nd intermediate-mass disky Es 
(boxes and filled circles respectively) are from iBender et alj jlffi)2|) and [Faber et al.l il997l). Cluste r Es (open circles) are 
taken from lPahrd Jl999l) and LIRGs ( diamondsj^om lshie^fcFischSl J 19981) and Tjamesetal] Jl999t) . The ULIRGs and the 
IR-bright QSOs of our program appear as triangles and filled stars respectively. The RL/RQC QSOs _of|Dunl2^^^^Jj2^0^) 
are open stars. Their black hole masses are converted into dispersions using the Mbh — G relation flVemain^etalT I^OO^)^ 
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Fig. 2. PG 1426+015- On the left panel the acquisition image with the two nuclei is shown. On the right panel, the NIR 
spectra are presented; the on-source integration time on the QSO is ~ 2 times deeper than that on the secondary nucleus. 
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